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ELECTRON  TRANSMISSION  THROUGH  AND  REFLECTION  BY  AN 


ARBITRARY  SET  OF  SQUARE  POTENTIAL  STEPS 


Dr.  Forrest  L.  Carter 
Chemistry  Division 
Code  6171 


INTRODUCTION 

Traditionally  solutions  of  the  Schrodinger  wave  equation  for 
square-well  and  step-function  potentials  have  served  primarily 
heuristic  and  pedagogical  purposes.  More  recently,  however, 
advances  in  the  use  of  multibeam  evaporators  have  made  possible 
the  growth  of  layered  semiconductors  with  layers  only  a  few 
Angstroms  thick.  In  addition  we  have  proposedfll  that  analogues 
to  NAND  and  NOR  switches  might  be  possible  at  the  molecular  size 
level.  That  proposal  relies  on  the  quasi-classical  work  of 
Pschenichnovf^ ]  who  for  a  one  dimensional  system  of  identical 
potential  barriers  anticipated  the  existence  of  an  electron 
tunnel  effect. 

The  tunnel  effect  occurs  when  the  electron's  energy,  co, 
matches  that  of  a  pseudo-stationary  state,  En,  (i.e.,  a  virtual 
level)  related  to  an  energy  level  of  the  well  between  barriers, 
when  a)  is  less  than  the  height  of  the  barriers.  The  work  of 
Pschenichnovf 2 ]  suggests  that  where  this  condition  is  met  (i.e., 

<d  =  En)  then  the  transmission  coefficient  T  is  equal  to  1;  in 
other  words,  at  this  energy,  the  barriers  are  perfectly  trans¬ 
parent.  However,  at  a  small  energy  mismatch,  transmission  falls 
exponentially  with  the  number  of  barriers.  The  extremely  rapid 
decline  in  electron  transmission  with  the  energy  mismatch 
provides  the  basis  for  a  switching  mechanism  assuming  that  one  is 
perturbing  one  or  more  of  the  potential  wells  or  barriers  and 
hence  disturbing  the  energy  level  of  the  pseudo-stationary  state. 
Both  the  layered  semiconductors  and  the  molecular  tunnelling  as 
conceived!!)  can  be  reasonably  treated  as  a  one  dimensional 
problem  (the  approach  taken  here) .  In  addition,  the  results 
obtained  below  may  be  useful  for  understanding  electron  transport 
and  scattering  along  a  polymer  chain  with  a  non-periodic 
distribution  of  substituents  or  conjugation. 
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Derivation 

From  the  author's! 3]  earlier  result  involving  a  periodic  set 
of  an  arbitrary  number  of  potential  step  functions  per  unit  cell 
it  seemed  reasonable  that  an  exact  expression  might  be  achievable 
for  electron  transmission  past  an  arbitrary  number  of  square  wells 
and  barriers.  That  possibility  is  vindicated  in  this  paper;  both 
the  derivation  and  the  generalized  analytic  expression  are  given 
for  the  transmission  coefficient  T  of  an  electron  past  any  number 
of  step-function  potential  changes.  The  method  is  classical 
textbook  (see  Schiff  [4]). 
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Fig.  1.  Square  well  (and  barrier)  potentials  of  width  d^  and 

depth  (height)  V£  in  the  path  of  a  unit  wave  approaching 
from  the  extreme  left. 


Before  developing  the  expression  for  perfect  transmission,  it 
is  very  useful  to  understand  the  relationships  for  the  wave 
function  in  the  regions  of  adjacent  potentials.  The  time 
independent  Schrodinger  wave  equation  in  atomic  units  is 


+  2[E  -  V(x) ]Y  =  0 


Eq.  1 


where  for  constant  potential  V  the  wave  function  ¥(x)  can  be 
readily  written  as 


fix)  =  A„e18nx  ♦  B  e-i6nx 
n  n  n 


Eq.  2 


correspondingly  to  waves  of  energy  E  *  oj  moving  to  the  right  and 
left  with  amplitudes  An  and  Bn  respectively.  The  exponential 
multiplier  Bn  is  determined  by  substituting  4^  (Eq.  2)  in  the 
differential  wave  equation  (Eq.  1) ;  while  the  coefficients  An, 

Bn  are  determined  by  joining  and  its  first  derivative  dV^/dx 


•  -i  m'*/'. 


continuously  at  each  step  position  a^  of  Fig.  1.  With  ¥n(x) 
corresponding  to  the  wave  function  with  an  £  x  <  an+i  the 
continuity  requirements  give  two  equations  for  each  potential  step 
at  the  positions  an;  thus, 

functional  continuity 


A  ,eiBn-lan  +  B  1e‘i0n-lan  =  A  e10nan  +  B  e_i6nan  Eq.  3 

n- 1  n-1  n  n 

slope  continuity 


iS  , a  i a  iB  a„ 

ft  a  o  n~l  n  Q  r,  -  n— 1  U  g  .  n  n 

3n-lAn-le  6n-lBn-le  PnAne 


Eq.  4 
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Obviously,  An  and  Bn  can  be  expressed  as  a  matrix  multiplication 
Mn  of  the  coefficients  An_i  and  Bn-i 


This  equation  can  be  easily  generalized,  as  was  done  by 
PschenichnovI 2]  for  his  quasi-classical  approach,  to  the  form: 


Fq.  6 


where  n  is  the  counting  index  of  the  potential  steps.  Clearly 
since  a  well  of  finite  barrier  requires  at  least  two  potential 
steps,  the  matrices  must  be  grouped  at  least  by  pairs  to  be 
compared  with  Pschenichnov' s  results. 

New  Results 


The  matrices  Mn  and  their  successive  products,  however,  are 
complex  (in  both  senses)  and  are  awkward  to  write  down.  It  is 
more  convenient  to  discuss  their  products  like 


n  n 

IIM.  and  n  M, 
1  1  k+1  1 


in  terms  of  the  new  and  much  simpler  analytic  expressions  U*Q 
and  V^D  to  be  given  explicitly  below  for  the  first  time  for 
various  n  up  to  6 . 

n  n 

In  these  terms  IIM.  and  H  M.  may  be  related  to  Eqs.  7  and  8 
1  1  k+1  1 

and  Eqs.  9  and  10,  respectively.  The  relationship  between  the 
coefficients  An,  Bn  and  A0#  Bft  in  terms  of  Ufiv  and  etc.  is 
indicated  by  Eqs.  Si  and  12.  °  *° 

The  and  VfjQ  are  to  be  given  in  Table  1  in  terms  various 
B^,  Si,  ana  Ci  as  defined  below  where  on  the  left  the  electron 
energy  oj  is  greater  than  the  potential  V^,  etc: 


Low  Potential 

High 

Potential 

V.  <  to 

1  — 

Vi 

> 

b) 

Bi  *  /ai  Vi 

*i 

✓V.  -  co 

l 

Ci  =  cos 

Ci 

■ 

cosh 

Si  »  sin  Bjdj^ 

Si 

- 

sinh  B^d^ 

Although  Ufj0  and  VpQ  functions  rapidly  become  rather  complex  for 
increasing  n  in  fact  tney  are  rather  simple  once  understood,  and 
they  can  readily  be  generalized.  The  verbal  algorithm  for  doing 
so  will  be  indicated  after  the  formula  for  the  transmission 
coefficient  T  is  developed. 


In  terms  of  Fig.  1  and  the  coefficient  pairs  Ac,  Bc  and 
An,  Bn  the  calculation  for  T  is  readily  described.  Approaching  the 
first  potential  step  at  Ai  from  infinity  is  the  unit  wave 
e*Bx,  thus  A0  *  1.  Reflected  from  that  first  potential  step  is 
the  wav  with  unknown  amplitude  B  *  B0.  The  coefficient  of  the 
wave  that  gets  through  to  the  region  V4,  x  >  a4  is  the  desired 
transmission  amplitude  T.  Mote  that  there  is  no  wave  coming 
from  the  right  in  this  region;  hence  An  ■  T  and  Bn  *  0.  From 
Eqs.  7  and  8  and  the  coefficients  just  specified  one  can  then 
readily  obtain  Eq.  13 


.‘"o'r'n'n1 
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no,2*IVno,2‘<0no,2"<Vno,2l  E|J-  13 
no  no  no  no  J 


U+  -iV+ 
no  no 


GENERATE  GENERAL  RULE  FOR  COEFFICIENTS 


TABLE  I  OF  COEFFICIENTS  IT  AND  V" 


TABLE  I  OF  COEFFICIENTS  IT  AND  V  (CONT’D) 

no  no 


BqB?Ba  ±  Bi  BiBs 


t 


then  with  A0  -  1  we  have  for  the  transmission  coefficient 


)2+(v+  ) 2-  (u“  )2-(v'  )2 

no  no  no  no 

4'Ko>2+,vno,2l 


1! 


Eq.  14 


Equation  14  informs  us  that  the  transmission  is  maximized 
when  u£0  and  v£0  approach  zero  while  the  numerator  remains 
finite.  Table  2  indicates  the  behavior  of  the  numerator  for 
the  first  few  values  of  n,  the  number  of  step  potential  changes. 
Generally  it  is  not  a  fast  varying  function. 


The  relationship  of  the  function  (obtained  here)  to  the 

solution  of  the  wave  equation  for  a  periodic  set  of  square  well 
potentials  is  rather  interesting.  This  is  readily  illustrated 
for  the  Konig-Penny  square  well  modelt3)  when  expressed  in  the 
notation  above.  Given  that  k  is  the  wave  vector  and  c  is  the 
unit  cell  length  we  have  (3J: 


cos 


kc  *  C^C 


1 

2 


S1S2 


Eq.  15 


This  may 


be  compared  readily  with  U*Q; 


(1  +  B3,C1C2~ 


Eq.  16 


/ 


Periodicity  after  two  potential  steps  requires  that  60=63  so  that 
one  may  clearly  combine  eqs.  15  and  16  to  form 

cos  kc  =  U^q/2  Eq.  17 

Given  n-1  steps  per  unit  cell,  Eq.  17  may  be  generalized  to 

cos  kc  =  U*  /2  Eq.  18 


Remembering  elementary  band  theory  we  note  that  allowed  bands 
exist  only  when  cos  kc  £  1  and  forbidden  gaps  occur  otherwise. 

Now  when  U^0  =  0  then  a  virtual  state  exists,  which  may  be  related 
via  Eq.  18  to  the  case  of  cos  kc  =  0,  i.e.,  the  center  of  the 
allowed  bands. 
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TABLE  2 

TRANSMISSION  COEFFICIENT  NUMERATOR 


N  =  (U+  )2  +  (V+  )2  -  (U-  )2  -  (V'  )2 
no  no  no  no  no 


N 


=  4  B0 
10  *  ^ 


N 


20 


«  %  <  ♦  si2» 


N30  "  4 


*[■ 


C  2C  2  +  s  2S  2  +  S  2 r  2 
C2  bl  s2  +  S2  C1 


+  SlV] 


This  then  agrees  well  with  one's  expectations  and  to  that  extent 
supports  the  conclusions  of  Pschenichnov  for  the  tunnelling  of 
electrons  through  periodic  barriers. 


Algorithm  for  U*_  and  Vr 
— •* - no - no 

+ 

Table  1  shows  that  in  U^q  only  even  numbers,  m,  of  S^  occur; 

further  that  the  number  of  such  terms  is  (n_1)  =  ■ —  that 

m  (n-l-m)ra! 


is,  all  possible  combinations  of  n-1  positions  take  m  at  a  time. 

B.  Bo  B5 

The  sign  in  front  of  a  term  like  (g—  +  C^C2C3S4S5  Is  +  1- 

5  p4  6 


mod  (m  DIV  2,2)  is  zero  and  negative  otherwise  (where  DIV  means 
integer  divide).  In  the  8^  ratios,  only  8i  corresponding  to  the 
Si  occur.  In  the  lead  8  ratio,  the  first  Ri  occurs  in  the 
numerator,  the  second  in  the  denominator,  and  the  alteration 
continues.  The  second  8i  ratio  is  the  reciprocal  of  the  first 
ratio  multiplied  by  80/®n’ 


Two  differences  occur  in  the  case  of  the  V^o  function. 

First,  m  is  only  odd.  The  second  occurs  in  the  $i  ratios.  In 
the  first  ratio  one  always  starts  with  80,  whereas  in  the  second 
ratio  the  first  possible  one  is  8i  and  the  last  8i  factor  is 
always  8n. 
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Summary 


Using  a  time  independent  method  the  analytic  formula  for  the 
transmission  of  an  electron  through  an  arbitrary  number  of  step 
potential  changes  has  been  derived.  An  algorithm  has  been  given 
for  writing  down  the  formula  for  any  number  of  potential  steps. 
In  addition  the  virtual  levels  have  been  identified  with  the 
allowed  band  centers  in  a  corresponding  case  of  n-1  potential 
step  changes  per  unit  cell  of  a  periodic  array. 

The  existence  of  a  generalized  analytic  solution,  now 
available  for  the  first  time  is  important  not  only  in  terms  of 
investigating  the  details  of  electron  tunnelling  in  one  and 
two  dimensional  molecular  switches  but  also  in  terms  of  electron 
conduction  in  polymers.  Thus,  while  eq.  18  is  applicable  to 
crystalline  or  periodic  polymers,  Eq.  14  will  be  useful  for 
understanding  and  modeling  non-crystalline  or  non-pcriodic 
materials  as  might  occur  in  multiple  co-polymers  in  the  presence 
or  absence  of  a  small  electric  field. 
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